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Abstract. This paper studies the structure of Lax pairs associated with integrable lattice 
systems (where space is a one-dimensional lattice, and time is continuous). It describes a 
procedure for generating examples of such systems, and emphasizes the features that are 
needed to obtain equations which are local on the spatial lattice. 
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I. Introduction. 



There has long been interest in integrable differential-difference equations (integrable lat- 
tice systems), especially since the discovery of the Toda lattice 1 . Such systems can have 
direct applications, for example in condensed-matter physics; and also occur as spatially- 
discrete versions of integrable partial differential equations 2 . Associated with each inte- 
grable lattice is a Lax pair, involving a matrix L which "steps" along the lattice, together 
with another matrix V which generates the time evolution. The purpose of this paper is 
to investigate the structure of this Lax pair, and how it affects the nature of the associated 
integrable systems. 

Throughout the paper, we work with functions (fi n (t) which depend on time t, and on 
an integer variable n. Such a function will be written simply as ip, its dependence on t and 
n being understood; then <p + denotes <p n +i(t), and denotes <^ n _i(£). The symbol A 
denotes the forward-difference operator, i.e. Atp = <p + — <p. For example, the Toda lattice 
equation (in first-order form) is 



We shall take the integer n to be unrestricted (i.e. the lattice is infinite). Our primary 
interest is in systems which are local, in the sense that the time-derivative of a variable 
at site n equals some expression involving the variables at sites n — 1, n and n + 1: i.e. 
nearest neighbours only. 

The point of view adopted here is that a lattice equation is integrable if it can be 
written as the consistency condition for a linear system (Lax pair) of a suitable type. This 
involves two q x q matrices L and V, the entries of which are functions of a "spectral 
parameter" A, as well as of t and n. In what follows, we shall, for the sake of simplicity, 
restrict to the case q = 2 (i.e. 2 x 2 matrices). The linear system is 



where ^ is a column 2- vector (depending on A, t and n). The consistency condition for 




dt 




(1) 




(2) 



(2) is 



^-L = V+L - LV . 



dt 



(3) 
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The crucial feature of (3) is that it specifies the evolution only of L, and not of V; so 
in order to get a meaningful equation, V has to be determined in terms of L. In the 
next section, we shall see how this happens. The subsequent sections illustrate how this 
structure can be used to generate integrable lattice systems. We shall see how known 
examples fit into this framework; and as a new example, construct a system which couples 
lattice versions of the Heisenberg ferromagnet and the derivative nonlinear Schrodinger 
equation. 

II. HOW L(A) DETERMINES V(X). 

In order to analyse the structure of L(X) and V(X), one needs to impose some requirements 
on the way that they depend on A. Let us assume that L, L _1 and V are rational functions 
of A, with poles at constant values of A (that is, the location of the poles does not depend 
on t or to). This is not the only possibility: for example, there is the well-known case of 
the lattice Landau-Lifshitz equation 3 , which involves elliptic functions of A. But we shall 
restrict to the rational case. 

By making a Mobius transformation on A, we may ensure that V^(A) is finite at A = oo, 
i.e. that its poles occur at finite values of A. Furthermore, since (3) is homogeneous in L, 
we have the freedom to multiply L by a scalar function of A (not depending on t or to). 
We can use this freedom to ensure that L is a (matrix) polynomial in A which is nonzero 
at each of the poles of V. Let p denote the degree of L as a polynomial in A. 

Equation (3) determines the evolution of each matrix coefficient in the polynomial 
L(A) = AX P + BX P ~ 1 + ... + £); so at this stage it is a set of coupled equations for q 2 (p + 1) 
functions (with q = 2 in what follows). As was emphasized above, the matrix V has to be 
determined in terms of L, since (3) does not specify its evolution: let us now examine how 
this happens. 

Assume for the time being that the poles of V are all simple. So V has the form 

N 

y(A) = ^yW (A _ Aa) -i + y(o) ) (4) 

a=l 

where V^°\ . . . , V^ N > are matrices independent of A. The general idea is that is 
determined by a choice of gauge, whereas each for 1 < a < N is determined by the 
residue of (3) at the pole A = A a . Note that equations (2) and (3) are invariant under the 
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gauge transformations 

L^A+LA" 1 , ( 5 ) 
r«AVA-' + (l A )A-', 

where A is a nonsingular 2x2 matrix depending on t and n (but not on A). A choice of 
gauge involves the following steps: 

(i) choose a form for L(A) (an algebraic condition on the entries in the matrix L), such 
that a necessary condition for this form to be preserved under gauge transformations 
(5) is that A + = A; 

(ii) then choose any which is consistent with (i) and the evolution equation (3). 

As an example to illustrate how this works, choose the coefficient of X p in L(X) (i.e. the 
leading term) to be the identity matrix. Then the remaining gauge freedom is (5), with 
A independent of n as required. And the leading term in (3) gives V± 0) = V<-°\ so any 
choice = V(°\t) then fixes the gauge completely; to obtain an autonomous system 
of equations, one chooses to be a constant matrix. This is the most straightforward 
choice of gauge; for other gauges, will depend on the functions appearing in L(A), 
and this dependence is in general nonlocal as we shall see below. 

Consider, next, the pole at A = A a . Clearly the residue of the right-hand side of (3) 
at this pole must vanish, i.e. 

V (") L («) _ L (a)y(a) = (g) 

where = L(X a ). The idea is that this constraint determines V^ a \ Since is a 
non-zero 2x2 matrix, there are two gauge-invariant possibilities: the rank of could 
be either 2 or 1. The next section deals with the rank 2 case; thereafter, we concentrate 
on the rank 1 case. 

III. Rank 2 Case: Nonlocal Systems. 

If is invertible, then (6) is a difference equation 

y{a) _ L (a)y(a) L (a)-1 ^ 

which determines Vn^ in terms of (say) L^\ : L^*} 2 , . . . and V"i^. In other words, 
is a nonlocal function of the entries in L^ a \ 



4 



To obtain a simple example which illustrates this case, take p = N = 1. Without loss 
of generality, we may set Ai = 0. Write L = AX + B, where B is invertible; and choose a 
gauge by specifying 

A-(° 1 

and = 0. So V(A) = V^/A, where is the solution of the difference equation 

Vl 1] = BV^B' 1 . (7) 

In general, the 2x2 matrix B contains four functions; let us effect a reduction to one 
function y n (t) by taking B to have the form 

In order for the reduction to be consistent, we need to have the form 

yd) _ f b 



-oo 



c 



(plus a multiple of the identity, which has no eventual effect). The resulting system of 
equations for y is 

A y = 6exp(y)-cexp(-y), (9) 
where Y n (t) is given by the non-local expression 

n-l 

Y n = y n + 2 ^ Vk- 

k= — ao 

This equation can be transformed to a form which looks local, by writing y n = (f) n — (p n +i- 
In terms of <j> n (t), eqn (9) becomes 

^(0+ - </>) = cexp(0 + + 4>) - 6exp(-0 + - 0). (10) 

If (for example) 6 = 0, then (10) is a differential-difference version of the Liouville equation 

4> tx = exp(20). (11) 

To get (11) from (10), we interpret c as the lattice spacing, put x = nc, and take the 
continuum limit c — > 0. Similarly, putting 6 = c in (10) gives a differential-difference 



version of the sine-Gordon equation 4 . But for these lattice Liouville and sine-Gordon 
systems, x and t (or n and t) do not represent space and time; rather they are characteristic 
(null) coordinates. In particular, one cannot specify arbitrary initial data at t = 0. In these 
lattice equations, one of the characteristic coordinates has become discrete, while the other 
(namely t) remains continuous. 

This example can be generalized in several directions, as follows. If one does not make 
the reduction (8), then one obtains an equation for the matrix B. Choosing the slightly 
different gauge A = I (the identity matrix) and = 0, leads to 

^-B-AV^=0 (12) 
dt v ' 

together with (7). Now (7) and (12) ensure that there exists a matrix R n (t) with 

B = l + R+R-\ V™ =-(±R)R- 1 . 

Lit/ 

Then (12) becomes 

±(R + R- 1 ) + A[(±R)R- 1 ]=0. (13) 
This is a differential-difference version of the principal chiral equation 

{R x R- l )t + (RtR~ 1 )x = (14) 

in which, as before, one of the characteristic coordinates has become discrete. 

This chiral equation generalizes, of course, to larger matrices (q > 2). Similarly, the 
Liouville and sine-Gordon examples generalize to differential-difference versions of other 
Toda field equations. 

Finally, it might be noted that there are difference-difference versions of the principal 
chiral 5 and Toda field 6 equations in which both characteristic coordinates (here x and 
t) become discrete. Another, very general, example of this type is the Hirota bilinear 
difference equation 7 . 

So systems with det ^ may be thought of as time-evolution equations which are 
nonlocal on the spatial lattice, or as equations where a characteristic coordinate (neither 
space nor time) has become discrete. To get local evolution equations, it is necessary for 
each to have rank 1. From now on, we shall concentrate on this case. 

IV. Rank 1 Case: Local Systems. 
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Given that has rank 1, the constraint (6) may be solved as follows. Write = 
L)?'L( a \ Assuming that trK^ is non-zero, the general solution of (6) is 

V(a) = t^R [f {a)Kia) + 9 {a) adj K^] , (15) 
where adj denotes the adjoint matrix of K^ a \ and where and are scalar 
functions, with f^' = f( a \ So the constraint (6) does not determine V^-\ . . . , 
uniquely: in particular one has the arbitrary functions g^. However, there is a further 
constraint, namely that the condition detZ/ Q ) = has to be preserved by the evolution 
(3). This gives equations on the g^ a \ which are precisely that they are constant on the 
lattice, just as the are: g^ = g^ . Then (15) can be re-written as 

y(«) = ° a ] . K^ + d^I, 

where and are functions of t only. It is clear that the term will not contribute 
in the evolution equations, and so only the remain; we may as well set d^ = 0, and 
take 

y(«) = C<X ( (16) 
trA» v ' 

At this stage, the could still be functions of time t; for simplicity, let us take them 
to be constants. One point to note about (16) is that is local: Vn^ is expressed in 
terms of L n and L n -i- 

So to obtain local evolution equations with 2x2 matrix Lax pairs, one first specifies 
the integer p (the degree of L(X)); the integer N appearing in (4) equals 2p, since the A a 
all have to be roots of detL(A). The matrices V^ a ' for 1 < a < 2p are given by (16), and 
involve the 2p constants c^ a ' . Finally, there is the choice of gauge, which determines V^ ^ . 
In general, turns out to be nonlocal, and special gauge choices are needed to ensure 
that it is local. 

One can relate all this to the r-matrix description (see, for example, reference 3; 
and also reference 8 which addresses the construction of an r-matrix from a given Lax 
pair). Suppose one has an L(X), a Poisson bracket and an r-matrix such that the Funda- 
mental Poisson Bracket Relations are satisfied. Suppose also that there exist Ai, . . . , Xn 
such that detL(A Q ) = for each a. Let r(A) be the trace of the monodromy matrix 
. . . L2(A)Li(A)Lo(A)L_i(A) . . . (which propagates from n = — oo to n = +oo). Then 

N 

H=Y,c {a) logr(A a ) 

a=l 



is a local Hamiltonian, and the corresponding Hamiltonian equations are just (3); the 
constants are the same as those appearing in (16). The problem from this point of 
view is to choose L(X), in a suitable gauge, such that a compatible r-matrix structure 
exists. 

More generally, one wants r(A) to be conserved in time, for all A — this then gives 
infinitely many conserved quantities. If one has a Lax pair (2) and boundary conditions 
which imply that V +00 = V-oo, then r(A) is indeed conserved. When V(X) dependes locally 
on the fields, then then it is easy to ensure that this condition is met; if, on the other hand, 
V is nonlocal, then conservation of r(A) is not guaranteed. This is one reason why locality 
is desirable, in the present context. 

If p = 1 (in other words, L(X) is linear in A), then detL(A) is a quadratic polynomial 
in A, the roots of which are Ai and A 2 . For the time being, let us assume that these roots 
are distinct; and by translating A set Ai = 1 and A 2 = — 1. It follows that L(A) has the 
form 

L(A) = ±(A + 1)LW-±(A-1)L( 2 ), (17) 

where l/ 1 ) and are 2x2 matrices each having zero determinant. So the entries in 
L( a ) involve six independent functions of t and n (in effect, the requirement that Ai and 
A 2 be constant has reduced the number of functions in L(A) from eight to six). The two 
satisfy evolution equations obtained by expanding (3): these are 

fL L («) = V (°) L {«) _ L («)y(0) _ i 3 ^ (18) 

where S = -L^V^ - +L^V (2 K The gauge choice reduces the number 

of functions by four (since A contains four entries), and we end up with a system involving 
two functions. A number of examples of this type are described in the following section. 

V. Some p = 1 Examples. 

This section exhibits some systems of the type described in the previous section. Such 
examples are simple to generate; but before doing so, we should ask when two lattice 
equations are to be regarded as being "the same" . More specifically, is there an appropriate 
equivalence relation on the set of all such systems? Certainly such an equivalence would 
include gauge transformations in which A was constant; and strictly-local redefinitions of 
the functions appearing in L (i.e. the new functions at lattice-site n depend on the old 



functions at site n only). However, it is customary to allow more general transformations 
than just these. A well-known case is that of the Toda lattice (1). If one replaces (p by 



r = Aif- = if — if- , (19) 

then (1) becomes 

A_ r = e ^-2e r + e r -; (20) 

and this is regarded as simply another form of the Toda lattice equation. 

But any equivalence relation which admitted (1) and (20) to the same class, would also 
have to allow the transformations <p i— > A k tp for all integers k (negative as well as positive). 
If one allows such highly-nonlocal transformations, then ones ends up with rather few 
equivalence classes; in fact, one might as well transform to action-angle variables, and say 
that the "only" integrable lattice is linear. Clearly this is inappropriate. 

The point of this argument is to conclude that there is no useful equivalence relation 
on integrable systems of the type that we are considering (unless we insist that (1) and 
(20) are to be regarded as distinct). This means that the task of listing such systems in a 
systematic way is not really well-defined. The best that one can do is to exhibit examples, 
and indicate how they are related to one another. 

Example (i). Choose a gauge such that — = 21. This is the gauge which was 
mentioned as an example in section II. As was remarked there, the gauge is then fixed 
completely by specifying some V^°\t). The simplest choice is to set 

y(o) = o. Note that 

L^ 1 ) and must have the form 

L (1) =I + M, L (2) = -I + M, 

where M is trace- free and detM = — 1. So we may write M = f • u, where ui, 02 and 
03 are the Pauli matrices, and f = f n (t) is a unit 3- vector. The dot denotes the usual 3- 
dimensional Euclidean scalar product (and A below will denote the vector product). The 
evolution equation for f, derived from (18), is then 

= A[(l + f_ • f)~ Vf_ + fit + i/f_ A f)] , (21) 

where /i = \{c^ — c^) and v = |i(c^^ + c^) are constants. Equation (21), then, is an 
integrable equation for the unit- vector function f = f n (t). If the parameter v is non-zero, it 
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can be set to unity by scaling t; so the system effectively depends on the single parameter 
/j,. The case \x = is the "Lattice Heisenberg Model" 3 ' 9 , so-called because it has the 
equation of the Heisenberg ferromagnet as a continuum limit. Indeed, if we set v = 2/h 2 , 
/j = jj/h, and let h — > 0, then (21) becomes 

f t = £fx + f Ai xx ; (22) 

the Heisenberg model corresponds to (i = 0. A slightly different choice of gauge, namely 
one in which is a nonzero constant matrix, yields a lattice nonlinear Schrodinger 
equation (different from the one in the next example). This is the lattice counterpart of 
the well-known gauge equivalence of the nonlinear Schrodinger and Heisenberg systems. 

Example (ii). Choose a gauge such that 

^=(1 °), £ (2 »=-(° ?V (23) 



x u J ' \0 1 / 

where u and v are functions of t and n. The remaining gauge freedom has A being a 
diagonal matrix function of t only. Substituting (23) into the evolution equation (18) 
determines V^°\ as 



T/(0) _ 1 ( —C^VU- C^V + C^V- \ 

v -2y c (2) u + c d) u _ - c m uv _ ) 



plus a diagonal matrix function of t, which by the residual gauge freedom may be set to 
zero. In addition, (18) gives the equations for u and v, namely 

A u = i(c (2) A« + c (1) A«_ + c {l) vuu- - c^vuu+), 

dt (24) 

— V = i(c (1) At) + C {2) AV- + C (2) UVV- - C (l, 1iTO + ). 

This is exactly the Ablowitz-Ladik system 2 ; their L-operator is slightly different from the 
one presented here (it is, in effect, quadratic rather than linear in A); but it is easily seen 
to be equivalent. In particular, if we choose c^ 1 ) = 2i = — c^ 2 ), and impose the (consistent) 
reduction v = ±k*, then (24) reduces to a lattice nonlinear Schrodinger equation 

i-j-it = u+ — 2w + it_ =F uu*(u + + U-) . 
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Example (Hi). Here we choose a gauge such that L^> is constant. Without loss of gener- 
ality, we may take l/ 1 ) to be 

iW = 0o)- (25) 
The most general form for the matrix L^> is 

L (2) = (UV UUj\ 

where u, v and w are functions of n and t. The evolution equation (18) for l/ 1 ) implies 
that must have the form 

v^) = - c —( Q uw \ c(2) ( vu - o\ r271 

2 v^- / 2(vu_ +«;)V W y V D J ' 

where AA = (i.e. A depends only on t). The residual gauge freedom, i.e. that which 
preserves (25), is (5) with 



s(n,t)/' 

this has to be used to determine A and D, and to eliminate one of the three functions u, 
v or w. In fact, the role of f(t) is simply to fix A: let us choose A = \c^ l \ The remaining 
freedom now is 

u i-> g-\ 

v i-> g+ v 7 (28) 



w\-*g 1 g+w. 



Equation (18) for gives equations for u, v and w, one form of which is 
^-logv = D+ + ±c {1) {wv_/v - uv) + ^c (2) " + + VW+ . , 



^log^ = A^-ic( 2 )^^V 
it V + w J 

) = A(-ic«™_ + ic^^— Y 
V 2 vu- +w J 



dt 
d 



(29) 



dt [VV 



d i I \ r\ 1(1")/ / \ 1 C2) ww + W - 

— log (mo) = — D — ±c K >(u + w + /u — uv) — ± c ( > ■ 



dt ~ 2 '2 W (t; U _+«;) 

(Any three of these equations implies the fourth.) We see from (28) that in order to 
remove the remaining gauge freedom, i.e. fix g (at least up to a function of t), there are 
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three possibilities. Namely, we can specify either w, or uw, or v as a function of the gauge- 
invariant combination uv. This in turn will determine D, and hence V^°\ on eliminating 
the relevant variable from (29). For example, specifying uw as a function of uv will give a 
local formula for D. Similarly, v can be specified as any function of uv. But if we impose 
w = F(uv), then we need F to be an exponential in order to get a local expression for D. 

This illustrates the way in which some choices of gauge lead to a nonlocal expression for 
y(o). 

As an example, let us take the gauge w = constant. Choose w = — 1, and write 
u = — e x , uv = y. Then (29) reduces to the system 



±y = i c (DA(y_e**-) + ±c< 2 >a( V ——), 

at z z \y — expAx-J 

A r (2) 



(30) 



dt ' 2 v ' 2(y-expAx_)' 

This is a version of the relativistic Toda lattice 10 M . 
Example (iv). The gauge choice 

where k is a constant, gives 



' V -/ce^" x 



(0 ) _ c^ 1 ) / 1 - ke y -+ Ax - -e x \ c^_( e x - y 
S - 1 -ke y -~ x - J 2 \ke~ x - 1 - ke~ y + Ax - 

and again leads to the relativistic Toda system 10 ' 14 ' 15 , this time in the form 
A y = _i c (i) A . A ( e w-+ Aa '-) - I c ( 2 )A;A(e-^+ A ^), 

Lit' 

A x = _i c (i) e tf(i + ke Ax ) + \c^e- y {l + ke Ax ~). 

CJ.6 



(31) 



The limit k — > gives the Toda system (1); it is worth examining this in more detail. In 
order to get (1), one may replace the variables x and y by <p and ij), where 

e x = -v^e^, e y = -1 + v^V, 

and set c (1) = c (2) = -l/v 7 ^. Then the k -> limit of (31) is indeed (1). But since 
— > oo in this limit, we need to re-interpret the associated Lax pair. The way to 
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get a well-behaved limit is to replace A by 2X/\ / k. The roots of detL(A) now occur at 
A = ±^v^, and so in the k — > limit they coincide. In fact, when k = we have 




and detL(A) has a double zero at A = (cf. ref 3). The corresponding expression for V^(A) 
is obtained by taking the k — > limit after first subtracting a constant multiple of the 
identity matrix: this yields 



VI. Some p = 2 Examples. 

In the p = 2 case, L(A) has the form L(A) = AX 2 + BX + C, where A, B and C are 2 x 2 
matrices; so to begin with, one has twelve functions of n and t. The requirement that the 
zeros A a of the quartic polynomial detL(A) be constant imposes four equations on these 
functions, and choice of gauge imposes a further four, so one is left with four independent 
functions. In other words, the generic system in this p = 2, q = 2 case is a system of 
coupled evolution equations for four lattice variables. 

Reductions of such systems, so that fewer functions are involved, are of course possible. 
One example that has been known for some time is a lattice version of the sine- Gordon 
equation in which space is discrete and time continuous (by contrast with the version 
mentioned in section III). Here the L-operator has the form 16 ' 3 



where (p and r\ are functions of n and t, h is a constant corresponding to the lattice 
spacing, and /(</?) = (1 + |/i 2 cos y?) 1 / 2 . The resulting integrable lattice has sine-Gordon in 
"laboratory coordinates" as a continuum limit: if we replace n by x = nh and let h — > 0, 
then (f satisfies ptt ~ <Pxx + simp = 0. 

In order to obtain another example, let us choose a different gauge, namely A = I, 
V'°) = 0. On B and C we impose the four constraints tr B = 0, detC = / constant, 
tr(BC) = 0, and trC + det B = -2k constant. It then follows that detL(A) = X 4 -2kX 2 +l 
has constant zeros. The matrices B and C now involve four independent functions, and 





Xf(ip)^ 
i/i(A 2 e-^/ 2 -e^/ 2 ) 
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their evolution is given by 

4s = AQ, 

d (32) 
— — C = i?_|_C — C 

(XL 

where Q = E a ^ (q) . # = ~ E a A" 1 !/^, and the are constructed as in (16). There 
are four parameters, namely the c^ a \ 

We can get an idea of what this system represents by looking at a continuum limit. 
To keep things simple, we assign particular values to the parameters, and this leads to the 
following continuum integrable system. 

Let B and C be 3-vectors, functions of x and £, satisfying the constraints B ■ C = 
and C ■ C = 1. Their time-evolution is given by 

B, = (CA B x ) x - {(C A C x • B)C} X + ±{(B • B)B} X , 

(33) 

C t = (C A C x ) x + (C A C x • B)B A C + ±(B • B)C X . 

This is integrable: it has a Lax pair of the form ty x = , \I/ t = , where U = 
— |i(BA _1 + CA~ 2 ) • a corresponds to L(A), and V = Efc=i ^4A _fc is a limiting version of 
the V(A) of the lattice system. 

The equations (33) have two obvious reductions: if B = 0, then we are left with 
the Heisenberg ferromagnet equation for C; while if C is a constant unit vector, then B 
satisfies the derivative nonlinear Schrodinger equation. So (33) may be viewed as a coupled 
Heisenberg-DNLS system; and (32) is a spatially-discrete version of this coupled system. 

VII. Concluding Remarks. 

It is clear that the examples given above provide only a very small sample of integrable 
lattice systems. One may envisage a classification involving the three integers q (the size 
of the matrices L and V), p (the degree of L(X)) and r (the maximum order of the poles of 
V(X)). But in view of the remarks at the beginning of section V, a complete classification 
would require a way of dealing with the problem of equivalence. 

We conclude with a remark on higher values of r. One can take a given L(A) (thereby 
fixing q and p), and allow r > 1: in other words, higher-order poles in V(X). This leads to 
hierarchies of lattice systems, of which the r = 1 cases are the first members. So hierarcies 
also fit naturally into this framework. 
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